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Abstract 

As a model, the Pais-Uhlenbeck fourth order oscillator with equation of motion 

is a quantum-mechanical prototype of a field theory containing both second and fourth order 
derivative terms. With its dynamical degrees of freedom obeying constraints due to the presence 
of higher order time derivatives, the model cannot be quantized canonically. We thus quantize it 
using the method of Dirac constraints to construct the correct quantum-mechanical Hamiltonian 
for the system, and find that the Hamiltonian diagonalizes in the positive and negative norm 
states that are characteristic of higher derivative field theories. However, we also find that the 
oscillator commutation relations become singular in the toi — > uj2 limit, a limit which corresponds 
to a prototype of a pure fourth order theory. Thus the particle content of the toi = 102 theory 
cannot be inferred from that of the tui 7^ uj2 theory; and in fact in the uJi — s- uj2 limit we find that 
all of the uji 7^ uj2 negative norm states move off shell, with the spectrum of asymptotic in and 
out states of the equal frequency theory being found to be completely devoid of states with either 
negative energy or negative norm. As a byproduct of our work we find a Pais-Uhlenbeck analog 
of the zero energy theorem of Boulware, Horowitz and Strominger, and show how in the equal 
frequency Pais-Uhlenbeck theory the theorem can be transformed into a positive energy theorem 
instead. 
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INTRODUCTION 



While attention in physics has by and large concentrated on second order equations of 
motion, nonetheless, from time to time there has also been some interest in higher derivative 
theories, with a typical such higher order equation of motion being the scalar field equation 
of motion 

id',-\/'M-W' + M')<Pix,t) = , (1) 

a wave equation which is based on both second and fourth order derivatives of the field. 
With the propagator associated with Eq. (P) being given by 
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^2 f.2 + 



(2) 



in momentum space, the expectation of canonical reasoning (see e.g. jlj] for a canonical 
study of theories of such second plus fourth order type) is that a quantization of the theory 
associated with Eq. ([T]) would, for 7^ 0, lead to a l/k"^ — 1/(A;^ + M^) dipole spectrum 
consisting of two species of particles, one possessing positive signature and the other negative 
or ghost signature. And indeed, part of the appeal of such propagators is that precisely 
because of this ghost signature, the propagator has much better behavior in the ultraviolet 
than a standard second order l/k"^ propagator, to thus enable this higher order theory to 
naturally address renormalization issues such as those associated with elementary particle 
self-energies or with quantum gravitational fluctuations. 

While a similar conclusion regarding the presence of ghosts might be anticipated to apply 
to pure fourth order theories as well (viz. the = limit of Eq. (0)),j2|] as we see from 
the form of Eq. (j21), the l/M^ prefactor multiplying the l/k"^ — l/{k'^ + M^) dipole term is 
singular, with it thus not being reliable to try to infer the structure of the = spectrum 
from that associated with the 7^ one. In fact, it is actually not even reliable to try to 
infer the spectrum associated with Eq. (0) via canonical reasoning at all, since even when 
0, the theory is a constrained one due to the presence of higher order time derivatives; 
and thus before drawing any conclusions at all, one must first quantize the theory in a 
way which fully takes these constraints into account, and only then identify the particle 
spectrum. To address this issue we shall thus effect a full Dirac constraint quantization 

of a quantum- mechanical prototype of Eq. something which as far as we know has 
not previously been carried out in the literature. Specifically, we shall study a restricted 
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version of Eq. (^J in which we speciahze to field configurations of the form (j){x, t) = q{t)e 
configurations in which Eq. (0) then reduces 



,ik-x 




(3) 



where 
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(4) 



with Eq. ^ itself reducing to the equal frequency uji = 002 when M = 0. As such, for 
general Ui 7^ 002 Eq. Q thus serves as a quantum-mechanical prototype of a field theory 
based on second order plus fourth order derivatives of the field, while becoming a prototype 
of a pure fourth order theory in the equal frequency limit. As well as being a model which 
we can actually solve analytically, the quantum-mechanical prototype encapsulates many 
of the general issues associated with the quantization of the full field theory, since it is the 
time derivatives of the fields rather then their spatial derivatives which are of relevance for 
constructing conjugates of the fields. With the equation of motion given in Eq. Q being 
derivable from the Pais-Uhlenbeck Lagrangian 



(7 is a constant), we shall thus quantize the theory starting from this Lagrangian. [5] This 
will allow us to establish that the properly constructed Ui 7^ UJ2 Hamiltonian does indeed 
possess eigenstates of normal and ghost signature (states which themselves are then explicitly 
constructed), while also allowing us to monitor the uji 002 limit, a limit, which like the 
above limit, will be found to be highly singular. 

CLASSICAL DIRAC HAMILTONIAN 

For a quantization of the theory associated with the Lagrangian of Eq. we would like 
to treat q and q as independent coordinates, but cannot immediately do so since if q is to 
be an independent coordinate, we could not then use dL/dq as the canonical conjugate of 
q. To obtain a form which would be appropriate for quantization we therefore introduce a 
new variable x{t) to replace q, and compensate for doing so by additionally introducing a 
Lagrange multiplier X(t) and a substitute Lagrangian 




(5) 




(6) 
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with the Dirac constraint method assuring us that the Hamiltonian which is to ultimately 
emerge will then be independent of the Lagrange multiplier. With the Lagrangian of Eq. 
(jni) possessing three coordinate variables, g, x and A, we must introduce the three canonical 
momenta, p^, Pg and px, conjugates which for the Lagrangian of Eq. (jU)) evaluate to 

To implement the Dirac method, as well as use these relations to construct the Legendre 
transform of the Lagrangian of Eq. (jH)), viz. 

Hl = Pxi + PqQ + P\\-L , (8) 

we also introduce two primary constraint functions for the canonical momenta which involve 
either the Lagrange multiplier or its conjugate, viz. 

(t)l=Pq-\ , (j)2=PX ■ (9) 

And then, rather than use H^, the Dirac prescription is to instead use 

Hi=Hl + Ui(j)i + (10) 

as Hamiltonian, where Hi is thus given by 

2 

Hi = 7r + ?(^? + ^2)^^ - 7;^W2<f + Ax + ui{pq - A) + U2PX ■ (H) 

For this theory we define generalized Poisson brackets of the form 

_ dA dB dA dB ^ dA dB dA dB ^ dA dB dA dB ^^^^ 
dx dpx dpx dx dq dpq dpq dq dX dp\ dpx dX ' 

so as to obtain the canonical 

{x,Px} = {q,Pq} = {X,px} = 1 . (13) 

Similarly, given the definition of Eq. (fT^ . we find that the Poisson brackets of the constraint 
functions with the Hamiltonian Hi are given by 

{(pl,Hi} = 7U;^Ci;^g - M2 + 0l{01,Ml} + 02{01,M2} , 

{4>2,Hi} = -X + Ui + 4>i{(f)2,Ui} + 4>2{(p2,U2} ■ (14) 
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Consequently, both of the {0i, Hi} and {02, Hi} Poisson brackets will vanish weakly (in the 
sense of Dirac) if we set 

Ui = X , U2 = 'jufulq . (15) 
On thus imposing these two conditions, Hi is then replaced by a new Hamiltonian 

2 

H2 = — + -{uj^ + uj^jx - -uJiUJ^q + PqX + -fUJiUJ^qpx ■ (16) 

With respect to this new Hamiltonian the constraint functions obey Poisson bracket relations 
of the form 

{01,i72} = 7^?^2PA , {02,i^2} = {PA,i^2} = O . (17) 

Hence, finally, if we now set p\ = Q (so as to enforce {0i,if2} = 0), the resulting algebra 
associated with the four-dimensional q,pg,x,Px sector of the theory will then (as befits a 
fourth order theory) be closed under commutation using Poisson brackets defined via 

_ dA dB dA dB ^ dA dB dA dB 
dx dpx dpx dx dq dpq dpq dq 

with the requisite classical Hamiltonian which we seek then being given by 
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TT Px . I 7/ 2 I 2\ 2 7 2 2 2 /in\ 

H = — + PqX + -{uj^ + uj^Jx - -(jJi(jj2q , (19) 



and with the requisite Poisson bracket relations which define the classical theory being found 
to be given by 

{x,px} =1 , {q,Pq} = 1 , 

Px 

{x, H} = — , {q,H} = x , {px, H} = -pq - j{cul + uj\)x , {p^, H} = -fulu^q . (20) 
7 

With such Poisson bracket relations the canonical equations of motion thus take the form 

x=— , q = x , Px = -Pq- li^^l + ^\)x , Pq = lujluh , (21) 
7 

to then enable us to both recover Eq. Q and make the identification x = q in the solution. 
Additionally, in solutions which obey these equations of motion the Legendre transform 
L = pxX + Pqq — if is found to reduce to Eq. (jSj) just as it should. The Hamiltonian H of 
Eq. is thus the correct one for the fourth order theory, and it thus is the one which is 
to be quantized. 
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CONNECTION WITH THE OSTROGRADSKI HAMILTONIAN 



In considering the theory associated with Eq. ()19p. it is important to distinguish between 
the general Hamiltonian H as defined by Eq. (jl9|) and the particular value i^sTAT that it 
takes in the stationary path in which the equations of motion of Eq. (j2H) are imposed, 
with the great virtue of the Dirac procedure being that it allows us to define a classical 
Hamiltonian H which takes a meaning (as the canonical generator used in Eq. ((201)) even 
for non-stationary field configurations (i.e. even for configurations for which p^i + PqQ — H 
does not reduce to the Lagrangian of Eq. (jSJ). Thus, for instance, it is the Hamiltonian 
of Eq. (fT^ which defines the appropriate phase space for the problem, so that the path 
integral for the theory is then uniquely given by /[rfg] [dp^ [dx] [(ipa;]exp[i / dt^p^x + Pqq — H)] 
as integrated over a complete set of classical paths associated with these four independent 
coordinates and momenta. 

As regards the stationary value that the classical i^sTAT takes when the equations of 
motion are imposed, viz. 

HsTAT = -q --{uj^+u^jq --u^u^q - > (22) 

we note that not only is this particular Hstat time independent (see e.g. Eqs. and 
fl26|) below), it is also recognized as being the Ostrogradski [3] generalized higher derivative 
Hamiltonian associated with the Lagrangian of Eq. (jSJ, viz. the Hamiltonian 

■ dL ..dL .d (dL\ ^ 

which Ostrogradski showed to be time independent in solutions which obey the generalized 
Euler- Lagrange equation 

3 , (24) 

while also being its associated time translation generator |8( . Moreover, not only is the classi- 
cal -ffsTAT time independent in the unequal frequency solution q{t) = aie~*'^^* + a2e~*'^2* + c.c. 
where it evaluates to 

ifsTAT(t^i 7^ UJ2) = 27(u;i - ujl){alaiujl - ala2Lol) , (25) 

it is even time independent in the equal frequency uji = uj2 = oj temporal runaway solution 
q{t) = cie"*"^* -|- C2te~*'^* + c.c. where it evaluates to 

ifsTAT('^i = ^^2) = 470;^ (2C2C2 + iujclc2 - iujc*2Ci) , (26) 
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dq 
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with a runaway in time not leading to a runaway in energy. (I.e. there is an appropriately 

defined energy which is time independent in the temporal runaway solution.) 

The value we have obtained for ifsTAT(t^i = ^2) in Eq- f|^(j|) recalls to mind the zero 

I 

energy theorem derived by Boulware, Horowitz and Strominger '9] for fourth order conformal 
gravity, another typical fourth order theory. Specifically, Boulware, Horowitz and Strominger 
constructed a classical energy for the fourth order gravity theory, and found that it would 
vanish identically if the gravitational field solutions were required to be asymptotically fiat. 
The analog condition for the equal frequency Pais-Uhlenbeck theory (a prototype of a pure 
fourth order field theory) would be to demand the absence of runaway solutions in time, 
and thus to set C2 = 0, a condition which would then yield for Hstat{^i = ^2) none other 
than the value zero. As we thus see, it is the restriction to an analog of asymptotic fiatness 
which leads to a zero value for the energy in the equal frequency Pais-Uhlenbeck case, with 
there being no c^Ci type term present in Hstat{^^i = ^2) even though Hstat{^i 7^ ^^2) 
contains both a^ai and 0302 type terms. Since the ci mode can only make a contribution 
to Hstat{'^i = ^^2) when C2 is non-zero, we can anticipate that in the quantization of the 
theory to be presented below the ci mode will not give rise to a propagating on-shell state, 
with the eigenspectrum of the equal frequency quantum Hamiltonian possessing not the 
number of energy eigenstates states associated with a two-dimensional harmonic oscillator 
(viz. the dipole structure exhibited in Eq. Q), but rather possessing only the number of 
energy eigenstates associated with a one-dimensional one. A further interesting feature of 
the form we have obtained for iJsTAT(^^i = ^^2) in Eq. ()2fi|l is that there is choice of sign 
for the parameter 7 for which the pure C2 contribution to the energy, viz. 87^7^0202, is then 
positive definite. On its own then, the runaway mode does not give rise to any disease 
such as a negative energy, with it yielding a positive one instead. There is thus no need to 
require the absence of runaway solutions in the equal frequency Pais-Uhlenbeck theory (by 
demanding an analog of asymptotic fiatness 1^) since in and of themselves they give rise 
to perfectly acceptable energies, and do not thus need to be avoided. Since the runaway 
solution classical energies are well-behaved, we can thus anticipate that the quantum energy 
eigenspectrum associated with the C2 mode sector will not possess any energy eigenstates 
with either negative energy or negative norm, even while the structure found for the unequal 
frequency Hstat{^i 7^ ^^2) in Eq. (^3)) indicates that there will be energy eigenstates of 
either negative energy or negative norm in the quantization of the unequal frequency case. 



7 



CONSTRUCTING THE UNEQUAL FREQUENCY FOCK SPACE 



With the Hamiltonian of Eq. ()19|) being defined for both stationary and non-stationary 
classical paths, a canonical quantization of the theory can readily be obtained by replacing 
{i times) the Poisson brackets of Eq. (j2(J|) by canonical equal time commutators. However, 
without reference to the explicit structure of the Hamiltonian itself, we note first that the 
identification (as suggested but not required by the equations of motion 11]) 

q(t) = aie-*"^^* + ase'^'^^* + H.c. , Pg{t) = i-fujiujjaie~''^'^ + i-fujfuj2a2e-''^^^ + H.c. , 
x{t) = -iujiaie-''^'' - icJaOae-^'^^* + H.c. , p^{t) = -^ula^e-'"^^^ - ^ujla^e'''^'-' + H.c. 

(27) 

then furnishes us with a Fock space representation of the quantum-mechanical commutation 
relations 

= = ^ , = \x,v^ = [q^Px] = [Px,Pq] = (28) 

at all times provided that 

[ai,a{] = -\ 2- , [a2,4] = ^ ^ , [ai,4] = , [01,03] = . (29) 

(Here and throughout both ui and uj2 are taken to be positive.) Then in this convenient 
Fock representation the quantum-mechanical Hamiltonian is found to take the form 

H = 2'j{ujI - ujl){ujlalai - 002(12(^2) + ^(^^i + ^^2) (30) 



with its associated commutators as inferred from Eq. ()20|1 then automatically being satisfied. 
With the quantity 'y{ujf — 002) being taken to be positive for definitiveness, we see that the 
[02, 02] commutator is negative definite, and with H being diagonal in the Oi, 02 occupation 
number basis, we see that the state defined by 

ai\n) = , a2\n) = (31) 



is its ground state, [l^ that the states 



\ + l) = [2^uJi{ujf-ujl)Y/^a\\n) , \-l) = [2^uj2{ujl-ujl)Y/^al\n) (32) 

are both positive energy eigenstates with respective energies ui and UJ2 above the ground 
state, that the state | + 1) has a norm equal to plus one, but that the state | — 1) has norm 
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minus one, a ghost state. Thus, as anticipated, the correct Hamiltonian for the unequal 
frequency theory can be diagonahzed in a basis of positive and negative norm states, with 
the relevant negative norm state wave functions being explicitly constructed via Eq. (j32p and 
its multi-particle generalizations. With the eigenstates of H labeling the asymptotic states 
associated with scattering in the presence of any interaction Lagrangian L/ which might be 
added on to the original Lagrangian L of the theory, the effect of L/ would be expected to 
induce transitions between asymptotic in and out states of opposite norm, with the unequal 
frequency theory then being non-unitary, jl^ However, even though the Hamiltonian of the 
unequal frequency theory does have ghost eigenstates, since the commutation relations given 
in Eq. (|29|) become singular in the equal frequency limit while both the Hamiltonian of Eq. 
(j3(jp and the normalized | + 1) and | — 1) states develop zeroes, we will have to exercise some 
caution in trying to discover exactly what happens to the ghost states when this limit is 
taken. 



CONSTRUCTING THE EQUAL FREQUENCY FOCK SPACE 

To explore the equal frequency limit it is convenient to introduce new Fock space variables 
according to 

1 / , 2huj\ 1 / , 2huj\ 

a = ai(l + ^)+a2(l-^) , h = ^{a,-a,) . (33) 



where 



^^iy^i±^ , ^ (^1 - (34) 

2 ' 2 ^ ' 



These variables are found to obey commutation relations of the form 

[a, a)] = \ , [a, 1)^] = ^ , [6, a)] = ^ , [6, h^] = v , [a, 6] = , (35) 

where 

As such, the introduction of the operators a, a\ b, is initially nothing more than a 
rewriting of the original Ol, o|, (32) '^2 

operators. Their utility derives from the fact that 
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when expressed in terms of them the coordinate q{t) of Eq. ()27|) gets rewritten as 



, 2ibuj . 
[a — ojcos et sm et 



+ H.C. , (37) 



and thus has a well defined e — > limit, viz. 

q(t, e = 0) = e-*'^*(a - 6 - 2ibujt) + H.c. . (38) 

Similarly, the Hamiltonian of Eq. (jHUj) which takes the form 

H = 8^uh\a^a - b^b) + 8juj\2b^b + a)b + 6^a) + uj (39) 

in the new variables, then limits to 

H{e = 0) = 87cj^(26^6 + a)b + fe^a) + , (40) 

while the following commutators of interest have limiting form 

[H{e = Q),a^]=uj{a^ + 2b'^) , [if(e = 0), a] = -cj(a + 26) , 
[H{e = 0), b^] = Lub^ , [H{e = 0),b] = -oob , 

[a + b,a^ + b^ = 2/t , [a-b,a^ - b^ = -2fi , [a + b,a) - b'^] = Q , (41) 

where now /i(e = 0) = /i = 1/(870;^). Thus even while the relation between the (a, b) 
and (ai, 02) sets of operators is explicitly singular in the limit, the limiting prescription 
obtained by using the a and b operators of Eq. nonetheless leads to the e = operator 
algebra given in Eqs. (j4Uj) and (j4H) in which there are no singular terms at all. However, 
this is not the only way to take the e — limit, as we could instead take the limit not of 
the operators of the unequal frequency theory, but rather of the states. There are thus two 
possible ways to construct the equal frequency theory - we can either start with the e = 
operator algebra and construct a new Fock space for it from scratch, or we can construct the 
equal frequency states as the e — > limit of the unequal frequency Fock space states. With 
these two prescriptions not being equivalent, we shall explore both of them, and discuss first 
the limit of the operator algebra. 

Taking the limit of the operator algebra 

For the operator algebra limit, we use the e — limit of the a, , b, W operators as 
the dynamical variables, with use of these operators then enabling us to construct an e = 
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theory which is well defined. For this theory the e = Fock space is then built on the Fock 
vacuum \Q) defined by a\Q) = b\Q) = 0, (we use this basis since, according to Eq. (f^ . 
states built on the unequal frequency Fock vacuum in which ai\Q) = 0, 02]^) = become 
undefined in the e = limit), with the state \Q) being an eigenstate of H with energy u. 
Moreover, in this e — limit we see that the [a, and [b, commutators both vanish, 
with, as we shall see below, there actually being two ways rather than one in which the Fock 
space states can implement this vanishing depending on how the e — > limit is actually 
taken. However, before constructing either of these two ways, we note that because the 
effect of the Hamiltonian is to shift by 2b^ in Eq. 1)411) . and because (unlike the unequal 
frequency case) neither of the operators {a^^ ±b^) which diagonalize the Fock space basis acts 
as a ladder operator for the Hamiltonian, il^] we can anticipate that the eigenspectrum of 
H{e = 0) will be quite different from the unequal frequency eigenspectrum. 

Because b^ does act as a ladder operator for H{e = 0), it is possible to construct a one 
particle energy eigenstate, viz. the state b^Q) with energy 2u!, but if we try to normalize it 
(our first way to realize the e = Fock space) the vanishing of the [b, 6^] commutator would 
entail that this state would have to have zero norm. On its own a zero norm state (unlike 
a negative norm state) does not lead to loss of probabilit y, b ut its very existence entails the 
existence of negative norm states elsewhere in the theory. [l5| Thus, with the state a^il) also 
having zero norm, we immediately construct the states 

states which obey 

(+|+) = 1 , (-|-) = -l , (+|-) = . (43) 

However, unlike the orthogonal positive and negative norm states | ± 1) of Eq. ()32|) which 
are eigenstates of the unequal frequency Hamiltonian, this time we find that neither of the 
I ±) states is an eigenstate of the e = Hamiltonian since 

H{e = 0)|±) = 2uj\±) + 4ujb^\Q) . (44) 

Thus even while the e = Fock space possesses negative norm ghost states, this time 
they can only exist off shell (where they can still regulate Feynman diagrams) but cannot 
materialize as on shell asymptotic in and out states. With a similar situation being found 
in the two particle sector 16| where only the state {b'^)'^\il) with energy 3a; is an energy 
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eigenstate, we see that the only states which are eigenstates of H{e = 0) are those of the 
form (b'^)"'\il), with all such states having zero norm and positive energy. 13] 

Now at first it is quite perplexing that the e = theory possesses far fewer energy 
eigenstates than the e 7^ theory, possessing only the number of eigenstates associated with a 
one-dimensional harmonic oscillator rather than a two-dimensional one. The reason for such 
an outcome derives from the fact that while the normal situation for square matrices is that 
the number of independent eigenvectors of a square matrix is the same as the dimensionality 
of the matrix, there are certain matrices, known as defective matrices, for which this is not 
in fact the case. A typical example of such a defective matrix is the non-Hermitian, Jordan 
block form, two-dimensional matrix 

M^(; 1) . (45) 

Specifically, while this matrix has two eigenvalues both of which are real (despite the lack 
of hermiticity) and equal to 1 (the trace of M is equal to 2 and its determinant is equal to 
1), solving the equation 

c kh:h:) 

leads only to g = 0, with there thus only being one eigenvector despite the two-fold degen- 
eracy of the eigenvalue, with the space on which the matrix M acts not being complete. 13] 
For the case of interest to us here, the action of the unequal frequency Hamiltonian of 
Eq. (IHUj) on the one-particle states a'^lQ), b^Q) yields 

2uj L -I 

Hb^n)= ^\e^a^\n) + {Auj^ -e^p\n)] . (47) 
2uj '- 

For this sector we can define a matrix 

1 /Au^ + e^ Auj^-e^\ 
"^''-^[ , 4.= -.0 

whose eigenvalues are given as 2ti;-|-e and 2uj — e (the trace of M(e) is and its determinant 
is Aui'^ — e^). For such eigenvalues, energy eigenvectors which obey 

i/]2a;±e) = (2t<;±e)]2c<;±e) (49) 



(46) 
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are then readily constructed as 

\2uj ± e) = [±eaT + {2uj ^ e)^^] \n) . (50) 

As we see, as long as e 7^ 0, the two one-particle sector eigenvectors of the Hamiltonian H 
are distinct. However when we let e go to zero, the two eigenvectors in Eq. (j^Uj) collapse 
onto a single eigenvector, viz. the vector &^|^^), the two eigenvalues collapse onto a common 
eigenvalue, viz. 2uj, and the matrix M (e = 0) of Eq. (jlSj) becomes {2uj times) the defective 
one given in Eq. (j45j) . Consequently, while the dimensionality of the a, a\ b, 6^ based Fock 
space does not change in the e — > limit (i.e. the dimensionality of the Fock space remains 
that of a two-dimensional harmonic oscillator), nonetheless, in the limit the Hamiltonian 
becomes defective, with H{e = 0) thus possessing far fewer eigenvectors than H{e 7^ 0), and 
with the dimensionality of the space of energy eigenvectors of H{e = 0) actually being the 
same as that of a one-dimensional rather than a two-dimensional harmonic oscillator. I1Q| As 
we thus see, the e — limit is indeed highly singular. Finally, with Eq. (jTfjl reducing to 

H{e = 0)a^n)= 2uj [a^n) + b^n) 

H{e = Op\Q)= 2ujb^\Q) (51) 

when e = 0, we confirm again that H{e = 0) only has one eigenvector, the zero norm 6''' 1^2) 
with eigenvalue 2u!. 

As well as these zero norm states, the e = Fock space also possesses harmonic oscillator 
coherent states as well. To construct these particular states it is convenient to introduce 
new operators 

a = a + b , P = a — b (52) 

which obey 

[a,a^]=2il , [p,P^] = -2ix , [a,P^]=0 , [p,a^]=0 , [a,P]=0 . (53) 
In terms of these operators the Hamiltonian of Eq. (j4(Jj) may be rewritten as 

H{e = 0) = ^ [2a^a - (3 - (3^ a] + u , (54) 
with and P'^ then being found to obey the relations 

[H{e = 0),a^]=uj{2a^ - (3^) , [H{e = 0), (3^] = ooa^ , 

[H{e = 0),a'^ - (3^=uj{a^ - (3^) , [H{e = 0),a - (3] = -uj{a - (3) . (55) 
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In terms of these operators we construct the coherent state 

\g) = e-'f''\n) (56) 

where the vacuum is again the one which a and b, and thus a and (3, annihilate. Us- 
ing the well-known relation e^"'"'^ = e'^e^e~^^'^^^'^ = e^e'^e~^^'^^^'^ which holds when 
[A, [A, B]] = [B, [A, B]] = 0, it can readily be shown that \g) has positive norm, 
viz. {g\g) = (f^le^^^^e^^e*^^^^""^""^-* = 1/e. Similarly, using the well-known relation 
e^Be-"^ = B +[A,B] + [A, [A,B]]/2 + [A, [A, [A,B]]\/6 + it can also be shown that 
application of the Hamiltonian to this state yields 

Hie = 0)\g) = ^(at2 _ ^t2 ^ 2a^P^)\g) , (57) 

with the coherent state thus not being an energy eigenstate. This same analysis generalizes 
to any other state of the form of an exponential of a string of creation operators acting on the 
vacuum, since the commutator of H{e = 0) with the string yields another string consisting 
also purely of creation operators, a string which thus commutes with the original string of 
creation operators. Thus whether we use states with a definite number of particles (such 
a^''"6''""|i7)) or states with an indefinite number of particles (such as \g)), we find that only the 
states are energy eigenstates, with all the e 7^ negative norm states having moved 
off shell in the limit. Consequently, in the e = theory there are no energy eigenstates with 
negative norm or negative energy. 

As well as this, another interesting aspect of the e = theory it is that even while it 
is a fourth order theory, the e = theory only possesses the number of observable states 
that would be found in an ordinary second order theory. Consequently, in a field-theoretic 
generalization of the equal frequency Pais-Uhlenbeck fourth order oscillator, we would an- 
ticipate that there would only be one observable particle in the theory rather than the two 
that would be associated with the field-theoretic analog of the unequal frequency case as 
exhibited schematically in Eq. Replacing a pure second order field theory by a pure 

fourth order one would thus not be expected to bring about any change in the number of 
observable states. 

Since the above equal frequency theory possesses zero norm states, as a quantum- 
mechanical theory it is thus somewhat unconventional. Even though the physical viabil- 
ity of the theory is not in question since the theory has been shown to possess no energy 
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eigenstates with negative norm, nonetheless, it is still of interest to try to recast the the- 
ory in an alternate form in which the theory has a more conventional look to it. Thus we 
seek to reinterpret the zero norm states as positive norm ones by using an unconventional 
scalar product with which to define the norm. To this end, we note that in constructing 
the eigenstates of a Hamiltonian, all that is needed is the introduction of a Hilbert space on 
which the Hamiltonian is to act to the right on a set of ket vectors \ip), with there being no 
need to introduce the dual vector bra states of those ket vectors for this purpose, with the 
structure of the energy eigenspectrum thus not being in any way sensitive to the structure 
of the dual space vectors. [20| The choice of definition of a scalar product is thus independent 
of the structure of the energy eigenspectrum, and while one ordinarily defines the dual bra 
vector simply as the conjugate of other choices are possible, and all of them provide 
legitimate formulations of quantum mechanics. 21] 

Thus for our purposes here, we note that with both of the equal frequency theory a^|Q) 
and states being zero norm states if their respective conjugates are canonically defined 
as (f2|a, in terms of the operators a and (3 which diagonalize the Fock algebra, we 

instead define the dual of the state [3'^\VL) to not be the standard but to be minus one 
times it (viz. —{VL\I3) instead, while at the same time we continue to keep the dual of to 
be the standard {VL\a. With such a dual state definition, the dual oil)^\VL) = {l/2){a^ 
is given by (l/2)(f2|(a + (3), viz. by {Q\a, with the overlap of b^\Q) with its own dual then 
being given by (l/4)(fi|(a + — = fi = 1/87^7^, a norm which, for positive 7, 

is then positive rather than zero. And similarly, the dual of a^\VL) = (l/2)(a^ + is 
given by (l/2)(r2|(Q; — j3), viz. by {^l\b, with the overlap of a^Q) with its own dual then 
being given by {l/4:){Q\{a — + = /i = l/87Co''^, a norm which is then also 

positive for the same choice of sign of 7. Moreover, the overlap of the dual of b^Q) with 
a^Q) is given by (l/4)(f2|(a + + an overlap which is zero, while the overlap 

of the dual of a^Q) with b^Q) is given by {l/A){Q\{a — — an overlap which is 

zero also. With respect to this definition of scalar product then, the states a^Q) and b'^\fl) 
form an orthonormal basis in which all states have positive norm. Thus because we do not 
need to specify the dual vectors in order to construct energy eigenkets, we have flexibility 
in defining what we mean by a norm and a scalar product, and can thus make the equal 
frequency theory zero norm states become orthonormal positive norm states instead, and 
can do so without in any way altering the fact that b^Q) is an energy eigenket while a^\Q) 
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is not. 

In addition to the above formal construction of dual vectors, there is also an explicit 
operational way to implement it. Specifically, we note that before we seek to define a new 
scalar product, the standard scalar products are given by (i7|aa'''|i7) = 0, (f2|a6'''|n) = fi, 
{Q\ba^Q) = fi, {Q\bb^Q) = 0. We now introduce an operator C whose action on Q;"^|r2) 
and is of the form Ca^il) — a'^\^), C/9^|0) = —P'^\il). In consequence, its action 

on a^Q) and b^Q) is of the form Ca^Q) — b^Q), Cb^Q) — a^|Q), an action which thus 
interchanges a^\Q) and b^Q) with each other. With respect to this operator C we thus 
define scalar products which evaluate to {il\aC a'^ = jl, (illaCb^Vl) = 0, {il\bC a'^ — 0, 
{Qj\bCb'^\Q) = p.. Thus with {Q\aC being defined as the dual of a'''|r2), and with {Q\bC being 
defined as the dual of 6^ 1^2), we thus generate an orthonormal basis. Now we saw above that 
the dual of a^fl) is given as {fl\b while the dual of b^fl) is given as {fl\a. In this construction 
then we thus identify the dual of a^\Q) as {Q\aC — and the dual of b^Q) as {Q\bC — 
{Q\a. Since the role of C is to interchange a^\Q) and we can thus conveniently describe 

the a^|n), b^Q) 2-space in the language of Pauh spinors. Thus on defining a^Q) — \ t), 
b^l^l) = I I), we see that the action of C on a'\il) and b^il) can be written as C| t) = | |), 
C| i) = I t). Similarly, its action on a^n) = (| t) + | i)) and = (| t) - | i)) can be 

written as C{\ t) + | i)) = (| t) + | i)), C{\ t) - | i)) = -(| t) - I i))- In the at|0), b^Q) 
2-space then the operator C can be represented by none other than the Pauli matrix ai , to 
thus establish that the needed C does indeed exist. Then given this definition of conjugate 
states, if we define wave functions via ipai^) — {x\a^Q,), ipb{x) — {x\b^Q) with respective 
conjugates ipl{x) = {Q\aC\x), ipli^) = {^\bC\x), we are able to obtain the conventional 
orthonormality relations for wave functions, viz. / dxip*{x)ipa{x) = J dx{Q\aC\x){x\a^\fl) = 
{fl\aCa^Q) = ft, J dxijjl{x)i/jb{x) = 0, / dxipl{x)tpa{x) = 0, / dxipl{x)ipb{x) = /t. 

As regards the Hamiltonian of the theory, we note that with respect to our definition of 
conjugate states, its matrix elements evaluate to {Q,\aCHa^\Q) — 1/470;^, {Q,\aCHb'^\Q) — 0, 
{Q,\bCHa)\Q) = 1/470;^, {Q\bCHb'^\Q) = 1/470;^. As we see, with respect to this definition 
of conjugates, the Hamiltonian remains defective, with its diagonal elements still being real. 
With there thus still only being one energy eigenstate, viz. the original with which to 

characterize asymptotic in and out S-matrix states, the theory remains unitary not because 
H has became Hermitian with respect to the new definition of conjugate states, but because 
it has remained defective. While we have thus shown that one could equivalently either 
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work with an equal frequency theory with zero norm states or recast the theory as one with 
positive norm states instead, to conclude this paper we now present an entirely different 
approach to the equal frequency theory, one that will lead to positive norm states without 
any need to introduce non-canonical dual vector states at all. 



Taking the limit of the Hilbert space 

While we have seen that it is possible to define an unconventional scalar product with 
respect to which all of the e = theory zero norm Fock states then have positive norm 
instead, it is also of interest to see if we could construct an alternate e limiting theory 
in which all energy eigenstates would have a positive norm even when the conventional 
definition of conjugate states is used. To this end we note that since we obtained zero 
norm states in the above by first going to the e = limit of the operator algebra before 
constructing the states, to avoid such zero norm states an alternate procedure would be 
to first construct a set of states with non-zero norm, and then take the e ^ limit while 
holding the norms of these states fixed. With the algebra of the a and h operators already 
being defined in Eq. (jH^j) even prior to taking the e — > limit, we can thus construct a basis 
for the e 7^ Fock space via states built on a vacuum which obeys a\VL) = b\Q) = rather 
than on one which obeys ai\Q) = 0, 02!^) = 0, and then explore the e — > limit of those 
particular states. We thus build normalized e 7^ states such as (the choice 7 < assures 
the positivity of A and u) the one particle 

|l,0) = ^|fi) , |0,1) = ^|1^) , (58) 



the two particle 



|2.0> = i^|n) . |l.l) = ^-,_^|n> . |0.2) = l^|n) . (59) 

and so on. 

Even while this particular basis is not a basis of eigenstates of the e 7^ Hamiltonian, 
for the e 7^ theory this particular basis is just as complete a basis as the one built out 
of the e 7^ eigenstates. However, unlike the e 7^ energy eigenstate basis, remarkably, 
and crucially for our purposes here as it will turn out, every single state constructed via the 
e 7^ and operators is found to have positive norm. It is not however an orthogonal 
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basis since overlaps such as (1,0|0,1) = /i/(Az/)^/^ are non-zero, overlaps which actually 
become singular in the e ^ limit. Since this basis is complete we can use it to define a 
particular limiting procedure in which the normalization of each of the particle states in Eqs. 
fl58|) and and their multi-particle generalization is held fixed while e is allowed to go to 
zero. Then in such a limit we find from Eq. ()58p that a^Q) and b^Q) both become null 
vectors. However, despite this, we cannot conclude that the creation operators annihilate 
the vacuum identically in this limit since matrix elements such as [a, 6'''] = /i do not 
vanish in the limit (/i(e = 0) = /t 7^ 0), with product operator actions such as a acting on 
being singular. Instead, the creation operators must be thought of as annihilating the 
vacuum weakly (i.e. in some but not all matrix elements), but not strongly as an operator 
identity. Now since neither of the states a^Q) or b^Q) survives in the limit, H{e = 0) now 
has no one particle eigenstates at all. With the two particle states {a^)'^\Q) and {b'^)^\Q) 
also becoming null in the limit, the only two particle state which is found to survive in this 
limit is the positive norm state a^6^|r2)//i (since the [a,b^ commutator does not vanish), 
and, quite remarkably, in this limit the state is also found to actually become an energy 
eigenstate, viz. 

H(e = 0)^|fi) = ^a^b^\n) + '^(b^yin) ^ ^a^b^\n) , (60) 

with the energy of this two-particle state being equal to Stu, which is nicely positive. With 
this analysis immediately generalizing to the higher multiparticle states as well, we see that 
the only states which then survive in the limit are states of the form (a^6^)'*|f2), positive norm 
states which also become positive energy eigenstates in the limit, with the observable sector 
of the theory thus being completely acceptable. We thus recognize two limiting procedures, 
first taking the limit of the algebra and then constructing the Fock space, or first constructing 
the Fock space and then taking the limit of its states, with this latter limiting procedure 
being an extremely delicate one in which the only states which survive as observable ones 



are composite. Thus, to conclude, we see that in the equal frequency limit the quantum 
theory based on the fourth order Lagrangian of Eq. (0) is in fact completely acceptable 
(in fact, technically, what we have shown is that the presence of energy eigenstates with 
negative norm in the unequal frequency theory is simply not a reliable indicator as to their 
possible presence in the equal frequency case), and it would thus be of interest to see the 
degree to which our results here might carry over to full fourth order field theories. 

18 



One of the authors (P.D.M) wishes to thank Drs. E. E. Flanagan and J. Javanainen for 
useful discussions. The work of P.D.M. has been supported in part by the Department of 
Energy under grant No. DE-FG02-92ER4071400. 



* lmannheim@uconiivm . uconn. edul 
|davidson@bgumail.bg u.ac.iIl 

[1] K. S. Stelle, Phys. Rev. D 16, 953 (1977); Gen. Relativ. Gravit. 9, 353 (1978). 

[2] With pure second order theories not possessing ghosts, and with second plus fourth order 
theories possessing them, it might be presumed that the ghosts arise from the pure fourth 
order theory itself rather than through an interplay between the second and fourth order 
theories. However, whether pure fourth order theories are to possess energy eigenstates with 
negative norm is something which has to be investigated in and of itself, and we are not aware 
of any demonstration in the literature that pure fourth order theories actually do possess 
negative norm energy eigenstates. 

[3] P. A. M. Dirac, Lectures on Quantum Mechanics, Belfer Graduate School of Science, Yeshiva 
University, New York (1964). 

[4] A. Pais and G. E. Uhlenbeck, Phys. Rev. 79, 145 (1950). 

[5] The need to properly take constraints into account was first noted by us in an earlier ver- 
sion of this paper (P. D. Mannheim and A. Davidson, Fourth order theories without ghosts, 
hep-th/0001115 , January 2000 (unpublished)) where the Dirac quantization was first pre- 
sented, and also by Hawking and Hertog (S. W. Hawking and T. Hertog, Phys. Rev. D 65, 
103515 (2002)) who quantized the Pais-Uhlenbeck oscillator by Feynman path integration 
techniques. 

[6] Indication of the singular nature of this limit is already seen at the classical level, with the 
equal frequency solution q{t) = cie"*"^* + C2te~*'^* + c.c. to Eq. @ (where uji = uj2 = uj) 
having a time behavior quite different from that of the unequal frequency solution where 
q(t) = aie-^"^!* + a2e-*'^2t + ^.c. 

[7] M. Ostrogradski, "Memoires sur les equations differentielles relatives au probleme des 
isoperimetres" , Mem. Ac. St. Petersbourg, VI 4, 385 (1850). 

[8] D. Bak, D. Cangemi and R. Jackiw, Phys. Rev. D 49, 5173 (1994). 

19 



[9] D. G. Boulware, G. T. Horowitz, and A. Strominger, Phys. Rev. Letts. 50, 1726 (1983). 
[10] While asymptotic flatness is natural to second order theories where the exterior solution to the 
second order Poisson equation V^0(r) = g(r) is given by the falling potential cj){r > R) = —a/r 
where a = dr'g{r')r'^, for the fourth order V^(f>{r) = g(r) the solution is instead given by 
(P. D. Mannheim and D. Kazanas, Gen. Relativ. Gravit. 26, 337 (1994)) the rising potential 
(j){r > R) = —fi/r + where (3 = dr'g{r')r'^/6, 7 = — dr'g{r')r''^ /2. Since there is 
no reason to expect dr'g{r')r''^ to vanish for an arbitrary source, even in the static sector 
asymptotic flatness would not appear to be an appropriate restriction to place on fourth order 
theories. 

[11] In general, equal time commutators such as [q{t),pg{t)] = i can be satisfied at all times by 
an [a, a^] = 1 commutator defined via q{t) = ae*-^*^*) /2^/^ + H.c, pq{t) = ia)e~^^'^*^ /2^/^ + H.c. 
with the function f{t) being arbitrary. The occupation number Fock space can thus be defined 
independent of the structure of and even independent of any interaction terms that might 
also be added on to the free particle H. Thus in general the dimensionality of the Fock space 
basis need not be the same as that of the eigenspectrum of and their basis states need not 
be in one to one correspondence. 

[12] Identifying as the annihilator of the Fock vacuum would eliminate negative norm states 
but would leave H without any lower bound on the ground state energy. 

[13] In their study of the Pais-Uhlenbeck oscillator. Hawking and Hertog actually argue that 
interactions do not in fact lead to any such transitions, thus making the wi 7^ 002 theory 
viable. 

[14] It is 6^ itself which acts as a ladder operator, though does not. 

[15] The unequal frequency theory states | + 1) it | — 1) both have zero norm. 

[16] The zero norm states (&^)^|r2) and (a^)^|r2) can be combined into positive and negative norm 

states while the state a)b^\^)/fl has norm plus one. 
[17] A further unusual property of these states is that the one particle matrix element {Q.\q{t)h'^ \^) 

is given by fie'^^'' rather than by the coefficient of the b field operator in Eq. H38|) . with there 

thus being a mismatch between the second quantized states and the first quantized wave 

functions. 

[18] While we ordinarily use Hermitian matrices in quantum mechanics since a Hermitian matrix 
necessarily has real eigenvalues, there is, nonetheless, no converse theorem which would require 



20 



the eigenvalues of a non-Hermitian matrix to be complex. Rather, hermiticity is only sufficient 
to secure real eigenvalues but not necessary, as the matrix M of Eq. (45) for instance directly 
demonstrates. Indeed, if we take a diagonal matrix Mi and add on to it a second matrix 
M2 whose diagonal entries are all zero and whose only non-zero entries are all located on 
only one of the two sides of the leading diagonal, then no matter what specific values these 
non-zero entries actually take, the secular equations for Mi and M3 = Mi + M2 will be 
absolutely identical, with M2 not contributing to [AI^ — XI\ = |Mi — A/| at all. As such, 
we can think of matrices such as M2 as being divisors of zero, so that just as the addition 
of zero to an ordinary number does not change its value, the addition of a divisor of zero 
such as M2 to a diagonal matrix does not change its eigenvalues. While there is no change 
in the eigenvalues, there can nonetheless still be a change in the eigenvectors, with there 
being a possible reduction in the number of eigenvectors in certain cases. It is thus possible 
to have a non-Hermitian A^-dimensional Hamiltonian with N real eigenvalues but less than A'^ 
eigenvectors. Apart from our interest here in the application of this phenomenon to the equal 
frequency fourth order oscillator, in passing we note that if applied to neutrinos, it could result 
in the non-observability of any right-handed or sterile neutrinos which might accompany the 
ordinary left-handed ones. 

[19] This reduction in the number of energy eigenstates in the e — > limit is a quantum-mechanical 
refiection of the fact that the equal frequency classical energy ffsTAT('^i = ^2) given in Eq. 

becomes completely independent of the coefficient ci of the classical solution q{t) = 
cie~*'^*-|-C2ie~*'^*-|-c.c. when the C2 coefficient is set equal to zero, with ci thus not representing 
a dynamical degree of freedom whose quantum analog could materialize as an on shell state. 

[20] The authors are indebted to Dr. A. Smilga for a very helpful comment in this regard, and 
also would like to thank him for informing them of his own recent work on ghost issues in 
fourth order theories - A. V. Smilga, Benign vs malicious ghosts in higher-derivative theories, 
|hep-th/040723l1 (2004). 

[21] For further discussion of such quantum- mechanical issues, see CM. Bender, D. C. Brody and 
H. F. Jones, Am. J. Phys. 71, 1095 (2003). 

[22] A model in which a field has a positive frequency part which annihilates the vacuum strongly 
and a negative frequency part which annihilates the vacuum weakly in a way such that only 
certain multiparticle states survive would appear to be a possible candidate mechanism for 



21 



quark confinement, with quarks themselves then only existing off shell. 
[23] Within such pure fourth order theories one of particular interest is the fourth order conformal 
gravity theory, a theory which has been found capable (P. D. Mannheim, Astrophys. J. 561, 
1 (2001)) of readily resolving the dark matter and dark energy problems which currently 
challenge the standard second order Newton-Einstein gravitational theory. The conformal 
theory is based on the imposition of the local Weyl conformal invariance gp,v{x) e'^"'^^^ g^iui^), 
and thus has a unique gravitational action of the form Iw = —dg J d'^x{—gY/'^C\fj_yt^C^^^'^ 
where C'^^^^'^ is the conformal Weyl tensor and ag is a dimensionless coupling constant. With 
the general conformal gravity rank two gravitational tensor —{—g)^^^^SIw/^gfiiy = 2agW^'^ 
reducing to = - W^WKp^/G in a g^^ = tj^^ + h^,^ linearization around 

flat spacetime (here K^"' = /i^'^ - r/^'^/i"„/4 and W = r]f"'d'^da - d^'d''), we find in the 
conformal gauge d^g^^ —g^'^ gupd^g'^P / 4: = (viz. a gauge condition which is left invariant under 
g^iv{x) e^"'^^'^ g^iu{x)) that the source- free region gravitational fluctuation wave equation then 
reduces to ag{dQ — V^)^i^^'^ = 0. With this equation of motion being decoupled in its tensor 
indices, we see that each tensor component precisely obeys none other than the M'^ = limit 
of Eq. (^. Since the conformal gravity theory is power counting renormalizable (due to the 
dimensionlessness of ag), it would thus be of interest to see if the structure we have found for 
the equal frequency Pais-Uhlenbeck theory might carry over to the conformal gravity theory, 
since it might then be possible to construct a fully renormalizable, fully unitary gravitational 
theory in four spacetime dimensions, one which despite its fourth order equation of motion, 
would nonetheless only possess one on-shell graviton and not two. 



22 



